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ABSTRACT 


The  paper  describes  one  type  of  ablative  thermal  protective 
structure  form  associated  with  returnable  satellites.  It  carries 
out  discussions  of  ablation  mechanisms  and  calculation  methods 
associated  with  FL  composite  materials  and  silicone  rubber.  In 
conjunction  with  this,  it  carries  out  comparisons  with  flight 
tests  and  gives  concluding  opinions . 


I .  INTRODUCTION 

Ablative  thermal  protection  structures  are  forms  which  are 
efficient  and  provide  relatively  high  heat  protection.  They  have 
the  greatest  adaptability  and  are  often  used  in  returnable 
satellites . 

Due  to  the  limited  nature  of  ground  aerodynamic  heat 
simulation  tests,  as  far  as  studies  of  calculation  methods  for 
heat  protection  associated  with  ablative  thermal  protection 
structures  are  concerned,  they  then  become  extremely  important. 

China,  in  the  late  1960's  and  early  1970's,  then  began 
studies  of  calculation  methods  for  the  ablative  thermal 
protection  of  spacecraft.  In  conjunction  with  this,  it  obtained 
calculation  results  which  marry  up  relatively  well  with  ground 
heat  simulation  tests. 


*  Numbers  in  margins  indicate  foreign  pagination. 
Commas  in  numbers  indicate  decimals. 
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As  far  as  the  engineering  calculations  associated  with 
ablative  thermal  protection  structures  are  concerned,  at  the 
present  time,  there  are  two  types  of  calculation  methods  which 
are  commonly  used.  One  type  is  designated  as  a  decomposition 
surface  method.  Calculation  formulae  can  be  seen  at  1  in  the 
Appendices.  It  is  recognized  that,  after  satellite  structure 
ablative  materials  have  been  heated,  thermal  decomposition  takes 
place  at  a  single  temperature.  This  is  also  to  say  that, 
ablative  materials  take  this  temperature  to  make  a  boundary, 
dividing  the  zone  of  the  carbon  layer  from  the  zone  of  the 
original  material.  In  this  way,  it  is  possible  to  keep  away  from 
the  complexities  of  pyrolytic  or  heat  decomposition  dynamics. 
Moreover,  this  makes  the  calculations  simpler.  With  regard  to 
the  mass  loss  rates  mp  associated  with  pyrolytic  gases,  by 
contrast,  one  uses  the  heat  conduction  difference  between  the  two 
sides  of  heat  decomposition  or  pyrolytic  surfaces  in  order  to 
make 

precise  determinations.  The  classic  case  associated  with  this 
calculation  method  is  seen  in  Reference  [ 1 ] .  Another  type  of 
calculation  method  is  designated  as  a  zone  decomposition  method. 
The  calculation  formulae  can  be  seen  at  2  in  the  Appendices. 

That  is,  recognizing  them  to  be  ablative  materials,  after  being 
heated  up,  pyrolysis  or  heat  decomposition  takes  place  within  one 
zone  or  region.  This  zone  or  region  is  defined  from  the  two 
precise  determinations  of  temperatures  for  the  top  and  bottom. 
This  type  of  calculation  method  is  relatively  complicated.  As 
far  as  different  types  of  ablative  materials  are  concerned,  it  is 
necessary  to  make  corresponding  simplifying  assumptions.  Only 
then  is  it  possible  to  realize  the  calculations.  Classic  cases 
can  be  seen  in  Reference  [2].  Speaking  in  terms  of  numerical 
value  calculations,  the  two  types  of  calculation  methods  are  both 
parabolic  type  partial  differential  equation  sets  associated  with 
the  solution  of  one  dimensional  dynamic  boundaries.  The 
treatments  of  this  type  of  dynamic  boundary  are  frequently  a  good 
deal  of  trouble  to  calculate.  In  particular,  as  far  as  the 
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amounts  of  sensitivity  to  these  influences  are  involved,  if  one 
selects  initial  thicknesses  for  carbon  layers,  makes  precise 
specifications  for  decomposition  surface  (zone)  temperatures,  gas 
mass  loss  rates  mp,  and  boundary  movement  speeds,  etc.,  and,  if 
the  values  selected  are  inappropriate,  then,  in  all  cases,  it 
will  influence  the  accuracy  of  the  calculation  processes. 
Attention  should  also  be  paid  to  rational  selections  for 
difference  form  choices,  and  the  physical  parameters  of 
materials.  This  will  also  give  rise  to  influences  which  cannot 
be  ignored  on  calculation  results. 

Due  to  the  universal  expansion  of  computers,  it  is  already 
possible  at  present  to  use  the  FORTRAN  language  in,  micro 
computers,  to  realize  design  calculations  for  ablative  thermal 
protection  structures.  /27 

Thermal  protection  design  work  associated  with  ablative 
thermal  protection  structures,  at  the  present  time,  are 
completely  capable  of  being  carried  out  by  thermal  design 
personnel  themselves.  At  the  same  time,  it  is  possible  to  make 
design  personnel  conveniently  carry  out  optimization  designs  for 
thermal  protection  structures.  It  is  also  possible  to  begin  the 
development  of  research  work  for  the  whole  series  of  such 
influences  as  thermal  and  physical  characteristics  of  materials 
on  design  results. 


II.  ANALYSIS  OF  ACTUAL  CASES 

The  exterior  shape  of  satellites  as  well  as  ablative 
structures  are  as  shown  in  Fig.l.  The  status  of  heat  received  by 
them  as  well  as  ablative  structures  is  as  shown  in  Table  1. 

As  is  shown  in  Table  1.,  the  nose  section  of  the  satellite 
and  the  skirt  section  opt  for  the  use  of  "FL  composite  materials 
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+  metal  shell"  structures.  Base  sections  opt  for  the  use  of 
"silicone  rubber  ablative  materials  +  metal  shell"  structures. 
Below,  we  carry  out  an  elucidation  of  the  ablative  thermal 
protective  mechanisms  for  these  two  types  of  structures 
respectively. 


Fig.l  Schematic  Diagram  of  Satellite  Exterior  Shape  as  well  as 
Ablative  Structures 

Key:  (1)  Direction  of  Thermal  Flow  (2)  Nose  Section  (3)  Skirt 
Section  (4)  Base  Section 


TABLE  1  THE  VARIOUS  SATELLITE  COMPONENT  STRUCTURES  AND  THERMAL 
FLOW  DISTRIBUTION 
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1.  FL  COMPOSITE  MATERIAL  +  METAL  SHELL  STRUCTURE 

This  type  of  structure  is  the  exterior  surface  of  a  metal 
shell  and  a  definite  thickness  of  FL  composite  material  combined 
together  into  one  body. 

(1)  FL  Composite  Materials  Ablative  Mechanisms  [3].  The 
ablative  mechanisms  of  FL  ablative  materials  are  capable  of  being 
divided  into  the  several  areas  below: 

(1)  Thermal  Breakdown  or  Pyrolysis  of  Materials.  As 
far  as  the  reception  of  heat  by  materials  is  concerned,  when  they 
reach  a  certain  temperature,  one  then  begins  thermal 
decomposition  or  pyrolysis.  Thermal  decomposition  or  pyrolysis 
temperatures  are  determined  on  the  basis  of  heat  mass  loss 
curves.  After  thermal  decomposition  or  pyrolysis  of  materials, 
the  gases  produced  will  escape  from  the  surface.  When  pyrolytic 
or  heat  decomposition  gases  pass  over  carbon  layers,  they  will 
take  along  with  them  a  portion  of  the  amount  of  heat.  At  the 
same  time,  they  will  also  produce  mass  induced  emission  effects. 
This  effect  is  then  capable  of  lowering  the  effects  of  convection 
adding  heat.  The  entire  thermal  decomposition  or  pyrolysis 
process  is  a  heat  absorbing  or  endothermic  reaction. 

(2)  Carbonization  of  Materials.  After  the  heat 
decomposition  or  pyrolysis  of  materials,  this  carbon  layer 
possesses  a  relatively  small  coefficient  of  heat  conduction  and  a 
relatively  high  radiation  coefficient.  Because  this  is  the  case, 
besides  giving  rise  to  insulation  effects,  it  is  also  capable  of 
radiating  to  the  outside  large  amounts  of  heat.  The  heat  flow 
radiated  to  the  outside  is  Qs=e~a  'Ta\  1 
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Because  this  is  the  case,  the  higher  the  temperatures  of 
carbon  layer  surfaces  are,  the  larger  are  then  the  amounts  of 
heat  radiated  to  the  outside.  This  process  is  also  endothermic. 

( 3 )  Combustion  of  Surface  Carbon  Layers .  When  carbon 
layer  surfaces  reach  a  certain  temperature,  the  carbon  will  burn. 

The  process  of  combustion  is  exothermic. 

(4)  Surface  Recession.  Due  to  the  combustion  of 
exterior  surface  carbon,  surfaces  will  recede.  If  the  carbon  on 
the  surface  is  unusually  loose  or  porous,  under  the  effects  of 
aerodynamic  forces,  one  will  also  have  the  occurrence  of 
mechanical  denuding  or  corrosion,  making  the  surface  recede. 

(5)  Conduction  of  Heat  Inward. 

(6)  Pile  Up  of  Carbon.  After  combustion,  in  carbon 
layers,  one  will  have  the  occurrence  of  a  phenomenon  of  carbon 
build  up.  This  process  will  increase  the  density  of  carbon. 

(2)  Simplification  of  Calculation  Models.  The  points  of 
difference  between  decomposition  surface  calculation  methods  and 
decomposition  zone  calculation  methods  exist  only  in  the  former's 
taking  decomposition  zones  and  seeing  them  as  surfaces.  Because 
this  is  the  case,  with  regard  to  the  simplification  of 
calculation  models,  the  two  are  capable  of  being  unified  for 
purposes  of  description  as  follows  (as  far  as  the  points  of 
difference  are  concerned,  they  will  be  pointed  out  respectively): 

/28 


(1)  One  Dimensional  Assumptions.  It  is  recognized 
that  the  whole  of  the  ablation  and  conduction  processes  occurs  in 
a  one  dimensional  space  in  all  cases. 


6 


(2)  Thermal  Equilibrium  Assumption.  It  is  recognized 
that  thermal  decomposition  or  pyrolytic  gases  maintain  local 
thermal  equilibrium  with  residual  carbon  layers.  Between  them, 
there  are  no  further  chemical  reactions. 

(3)  Induced  Emission  Effect  Assumption.  It  is 
assumed  that,  as  far  as  the  mutual  interactions  between  thermal 
decomposition  or  pyrolytic  gases  and  boundary  layer  air  are 
concerned,  it  is  possible  to  use  air-air  induced  emission  effects 
in  order  to  make  descriptions  and  that  there  are  no  chemical 
reactions  between  them. 

(4)  Surface  Recession  Assumption.  It  is  assumed  that 
carbon  layer  surfaces  have  only  surface  recession  given  rise  to 
by  ablation.  It  is  not  necessary  to  consider  recession  which  is 
given  rise  to  by  other  factors  (for  example,  mechanical  denuding 
or  corrosion) . 

(5)  Decomposition  Zone  Calculation  Method  Assumption. 
As  a  prerequisite  for  the  assumptions  about  the  various  items 
above,  it  is  recognized  that  thermal  decomposition  or  pyrolysis 
is  carried  out  within  a  thermal  decomposition  or  pyrolytic  zone. 
The  thermal  decomposition  or  pyrolytic  zone  is  precisely 
specified  by  the  use  of  precise  determinations  of  the  upper  and 
lower  temperatures  of  the  zone  (see  Fig.2).  Within  thermal 
decomposition  or  pyrolytic  zones,  the  thermal  and  physical 
characteristic  parameters  associated  with  materials  can  be 
recognized  to  be  functions  of  temperature. 

(6)  Decomposition  Surface  Calculation  Method 
Assumption.  Given  the  prerequisites  of  the  four  assumptions 
((1)),  ((2)),  ((3)),  and  ((4)),  it  is  recognized  that,  between 
carbon  layers  and  the  original  materials,  there  exists  an 
interior  surface  defined  by  a  fixed  temperature.  In  this  way,  a 
surface  is  designated  as  thermal  decomposition  or  pyrolysis 
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surface.  Using  the  thermal  decomposition  or  pyrolysis  surface 
for  substitution  with  the  thermal  decomposition  or  pyrolysis 
zone,  in  this  way,  an  assumption  is  designated  as  the 
decomposition  or  pyrolysis  surface  calculation  method  assumption. 


(3)  Discussion  Relating  to  Decomposition  Surface 
Calculation  Methods.  When  the  thickness  of  FL  composite 
materials  is  25mm,  a  comparison  between  temperature  calculation 
results  and  satellite  flight  test  results  is  seen  in  Fig. 3. 


Fig. 2  Schematic  Diagram  of 
the  Decomposition  Zone 
Calculation  Method  Model 

Key:  (1)  Original  Surface 

(2)  Ablation  Zone 

(3)  Carbon  Layer  Zone 

(4)  Decomposition  Zone 

(5)  Original  Material  Zone 

(6)  Back  Wall  Structure 


Fig. 3  A  Comparison  of  the  Telemetry  Results  and  Calculated 
Results  for  a  Nose  Section  at  a  4>  =21°  Location 


Key:  (1)  Temperature  (2)  Flight  Test  Values  (3)  Design 
Calculation  Values  (4)  Flight  Time  Axis  (5)  Design  Time  Axis 
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Fig. 4  A  Comparison  of  Nose  Section  Carbon  Layer  Thickness  Design 
Results  and  Telemetry  Data 

Key:  (1)  Thickness  (2)  Telemetry  Results  (3)  Calculation 
Results 


Speaking  in  regard  to  nose  sections,  due  to  the  fact  that, 
during  the  processes  of  return,  thermal  flow  is  relatively  high, 
material  ablation  occupies  a  leading  role.  Using  decomposition 
surface  calculation  methods,  it  is  possible  to  obtain  results 
which  marry  up  relatively  well  with  flight  test  results  (the 
higher  thermal  flows  are  in  sections,  the  better  they  match  up). 
See  Fig. 4.  However,  from  Fig. 3,  it  is  possible  to  see  that,  as 
far  as  comparisons  of  flight  tests  and  calculation  results  are 
concerned,  calculation  results  for  temperature  tend  to  the  low 
side  on  exterior  surfaces.  They  tend  to  the  high  side  on 
interior  surfaces.  The  main  reason  for  the  calculation  results 
on  exterior  surfaces  tending  to  be  low  is  that  carbon  layer 
thermal  conduction  coefficients  used  in  calculations  are  greater 
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than  actual  thermal  conductance  coefficients.  Moreover,  the 
principal  cause  for  inside  surface  temperatures  tending  to  be 
high  is  one  created  by  the  use  of  decomposition  surfaces  to 
substitute  for  decomposition  zones. 

As  far  as  skirt  sections  are  concerned,  due  to  the  fact  that 
thermal  flows  are  relatively  small,  insulation  plays  the  key 
role.  In  this  type  of  situation,  the  effects  of  decomposition 
zones  can  then  be  ignored. 

(4)  Discussion  Regarding  Decomposition  Zone  Calculation 
Methods.  Calculation  models  associated  with  decomposition  zone 
calculation  methods  are  as  shown  in  Fig. 2.  Thermal  protective 
systems  as  a  whole  are  composed  of  a  layer  of  heat  protective 
material  and  several  types  of  different  back  wall  materials. 
Between  the  various  layers  of  materials  of  the  back  walls,  it  is 
possible  to  have  air  spaces.  It  is  also  possible  not  to  have  air 
spaces.  The  thermal  and  physical  parameters  associated  with 
various  types  of  materials  can  be  understood  to  be  functions  of 
temperature.  Other  assumptions  are  as  described  in  (2).  On  the 
basis  of  these  fundamental  calculation  formulae,  decomposition 
zone  calculation  method  programs  were  written. 

Using  the  programs  written,  calculations  were  carried  out 
for  nose  sections  and  skirt  sections.  The  results  of  temperature 
calculations  matched  up  when  compared  with  flight  test  results. 
See  Fig. 5  and  Fig. 6. 


2.  SILICONE  RUBBER  ABLATION  MATERIAL  +  METAL  SHELL 
STRUCTURES 

Silicone  rubber  ablation  coating  takes  silicone  rubber  as  a 
base.  It  adds  in,  at  given  ratios,  polyethylene  and  composite 
short  fibers  as  well  as  the  corresponding  binders  and  promoting 
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agents  for  its  formation.  Taking  the  coating  and  spreading  it  on 
the  exterior  surface  of  metal  shells,  one,  then,  forms  ablation 
coating  type  thermal  protection  structures.  This  type  of  /30 
ablation  coating  possesses  relatively  good  ablation  properties 
and  insulation  properties.  It  is  appropriate  for  use  in 


Key:  (1)  Temperature  (2)  Surface  Temperature  Flight  Test  Values 
(3)  Surface  Temperature  Design  Calculation  Values  (4)  Flight 
Test  Values  for  Location  9mm  Away  From  Original  Surfaces  (5) 
Design  Calculation  Values  (6)  Time 


Fig. 6  Skirt  Section  Surface  and  Back  Surface  Temperature  Curves 

Key:  (1)  Temperature  (2)  Surface  Temperature  Flight  Test 
Values  (3)  Surface  Temperature  Design  Calculation  Values  (4) 
Back  Surface  Temperature  Flight  Test  Values  (5)  Back  Surface 
Temperature  Design  Calculation  Values  (6)  Time 
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(1)  Silicone  Rubber  Ablation  Coating  Ablation  Mechanisms. 
Silicone  rubber  ablation  coatings  also  belong  to  the 
carbonization  class  of  ablation  materials.  Because  this  is  the 
case,  there  are  no  differences  in  principle  between  their 
ablation  mechanisms  and  the  general  run  of  carbonization  ablation 
materials.  However,  silicone  rubber  coatings  also  have  a  number 
of  special  points.  First  of  all,  during  ablation  processes, 
materials  will  show  the  occurrence  of  expansion.  This  is  shown 
clearly  by  tests.  When  the  coating  thickness  is  3-4  mm, 
generally  speaking,  there  will  be  expansion  of  2-4  mm  in  all 
cases.  Secondly,  carbon  layer  structures  are  also  relatively 
complicated.  When  heat  flows  are  relatively  large,  tests  clearly 
show  that,  on  surfaces,  there  will  be  produced  a  hard  shell  which 
has  a  loose  and  brittle  texture.  The  lower  surface  is  a  layer 
possessing  a  number  of  powdery  half  air  pocket  structures. 

Again,  lower  surfaces  are  expansion  layers.  On  the  lower  surface 
of  the  expansion  layers,  by  contrast,  are  intact  silicone  rubber 
coating  layers.  Looking  at  everything  as  a  whole,  the  portion  of 
structures  which  belong  to  carbon  layers  are  relatively 
complicated. 

(2)  The  Simplification  of  Calculation  Models.  Due  to  the 
fact  that  the  mechanisms  of  heat  exchange  associated  with  the 
half  air  pocket  structure,  which  is  formed  by  the  double 
splitting  of  silicone  rubber,  are  still  unclear,  as  a  result  of 
this,  calculation  assumptions  are  also  very  crude. 

Opting  for  the  use  of  decomposition  surface  calculation 
methods,  simplification  assumptions  are  the  same  as  the  various 
items  ((1)),  ((2)),  and  ((3))  in  (2)  discussed  above.  Besides 
this,  in  a  general  model  of  carbonizing  ablation,  surfaces  recede 
in  all  cases.  However,  in  the  processes  of  the  ablation  of 
silicone  rubber,  expansion  is  a  primary  factor.  Because  this  is 
the  case,  one  need  not  consider  surface  recession  which  is  given 
rise  to  by  the  combustion  of  carbon. 
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As  far  as  the  expansion  of  silicone  rubber  is  concerned,  one 
need  only  give  consideration  to  its  surface  temperature  functions 
(see  Fig.7).  For  calculation  model  see  Fig. 8. 


Fig. 7  Change  Curves  for  Expansion  Layer  Thickness  6  Following 
Along  with  Surface  Temperature  Tw 


m  (D 

(X) 


Fig. 8  Silicone  Rubber  Coating  Ablation  Calculation  Model 

Key:  (1)  Carbon  Layer  (2)  Original  Surface  (3)  Decomposition 
Surface  (4)  Original  Coating  (5)  Metal  Shell 
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(3)  Discussion  Concerning  Silicone  Rubber  Ablation  Coating 
Calculation  Methods.  The  section  after  decomposition  of  silicone 
rubber  is  generally  known  as  the  carbon  layer.  The  heat 
conduction  coefficient  for  this  section  is  very  difficult  to 
actually  measure.  Because  this  is  the  case,  the  key  to 
calculations  rests  in  how  the  thermal  conduction  coefficients  for 
carbon  layer  zones  are  selected.  On  the  basis  of  analysis  of 
results  from  surface  tests,  one  makes  the  assumptions  that 
follow:  that  is,  taking  carbon  layers  and  dividing  them  into  two 

sections,  inside  the  expansion  layer  close  to  the  decomposition 
surface,  it  is  filled  with  thermal  decomposition  or  pyrolitic 
gases.  The  coefficients  of  thermal  conductance  are  relatively 
low.  It  is  recognized  that  the  thermal  conductance  coefficients 
are  relatively  close  to  the  thermal  conductance  coefficients 
associated  with  air  in  the  same  range  of  temperatures.  This 
value  is  selected  as  0 . 84xlO-4kW/m» K.  The  upper  surface  of  this 
layer  is  a  loose,  porous  layer  of  carbon.  The  heat  conductance 
coefficient  will  be  larger  than  the  heat  conductance  coefficient 
of  the  original  material.  Values  are  chosen  as  1.5  times  the 
original  material. 

After  going  through  these  kinds  of  assumptions,  option  is 
made  for  the  use  of  decomposition  surface  calculation  methods. 

The  calculation  results  are  shown  in  Table  2 . 

The  metal  shell  back  surface  temperatures  measured  in  flight 
tests  are  TB  ^  470  K.  The  maximum  temperature  for  TB  calculated 
results  is  505  K.  Giving  consideration  to  the  influences  of  the 
heat  flows  sustained  by  the  base  sections  and  errors  in  the  set 
up  of  the  flows  as  well  as  other  factors,  it  is  possible  to 
recognize  that  calculation  results  and  measured  results  marry  up 
relatively  well. 
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TABLE  2  TEMPERATURE  DISTRIBUTIONS  AT  DIFFERENT  INSTANTS  FOR 
ABLATION  COATING  STRUCTURES 


(£)  fit  I  a]  t  ( s ) 

(K)  j 

ZfeWMW&B;  Ta  (K) 

to  =  159 

300 

300 

300 

tqmax  —  23  1 

775 

471 

300 

^Tnu*  =  233 

781 

474 

300 

t T Bmax  ~  268 

465 

505 

300 

^Tcm»x  =  289 

453 

495 

300 

(S)  ft.  Steals])  —  &  ft  £  ffi  S  $  Bl  [B]  1  <TBm, 

*  —  £  Fi  fY  ®  ft  ft  &  ft  fu] , 

Key:  (1)  Time  (2)  Surface  Temperature  (3)  Metal  Shell  Back 
Surface  Temperature  (4)  Insulation  Material  Back  Surface 
Temperature  (5)  Note:  tp  is  time  of  origin;  tg  max  is  maximum 
thermal  flow  time;  tT  max  is  maximum  surface  temperature  time; 
tTB  max  is  metal  shell  back  surface  maximum  temperature  time; 
trfc  max  i-s  insulation  material  back  surface  maximum  temperature 
time. 


III.  APPENDICES 

1.  DECOMPOSITION  SURFACE  ALGORITHM  CALCULATION  PROCESSES 
Coordinate  relationships  (see  Fig. 9): 


Fig. 9  Decomposition  Surface 
Algorithm  Coordinate 
Relationships 

Key:  (1)  Ablation  Layer 

( 2 )  Carbon  Layer 

(3)  Original  Material  Layer 

(4)  Insulation  Material  Layer 


(D 

■ 

X 

0  MM 

X 

@)mmm 

~ 

X' 

( tymmai 

X • 

>  * 
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Equations : 


X'=: 


x=[‘^dt 

Jo  P 

(1) 

J  op  —p  J  „  p 

(2) 

J  Op'  -p 

(3) 

("  =  X\ 

(4) 

a(^)— ' ■%- 


CX^Y^X+X) 

CX  YX^Y^X  +  X+X') 

JL(ir''dT'  V  »  „  dT"  — 

8y\  dY  )~p ‘Cp  ‘~dT~  u+i+r<Ki0+i'0+r0) 


37 


(5) 

(6) 

(7) 


Original  conditions : 


7(7,0)  =4(7)=^  (8) 


Boundary  conditions: 


Carbon  layer  surface: 


7 =X 

y4>+qK-e'0'TwA+mc'Mc+k^-  = 


37 


Thermal  decomposition  or  pyrolytic  surface: 

Y-  X  +  X 

cT  =  T, 


When  Tr--i+x<Tt 


(9) 


(10) 
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m,=  0 

The  common  boundary  surface  between  original  materials  and 
insulation  materials  is  Y  =  X  +  X  +  X1 


(ID 
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ypf  J1// 

_k,dT'  _  „dT" 
d Y  k  ~§Y~ 


(12) 


The  insulation  material  back  surface  Y  =-.x  +  X  +  X'+  X'  ' 


dr" 

37 


=  0 


or 


oT" 

-k"-W=£,/'° •  C -r.o 


(13) 


2.  DECOMPOSITION  ZONE  ALGORITHM  CALCULATION  PROCESSES 


(see  appended  Fig. 10) 


Coordinate  relationships 


CO  ftse 

X 

(5)  m 

X 

(S') 

X' 

X’ 

X" 

\y 


X=X0+\'?±dt 

J  0  P 

(14) 

-  - 

(15) 

X  —  Y  t-t  g  i  —  X 

(16) 

X'  =Yt.ts1-  x-x 

(17) 

X"=X\-  ( X+X+X ') 

(18) 

X"’=X"t 


Fig. 10  Decomposition  Zone  Algorithm  Coordinate  Relationships 


Key:(l)  Ablation  Layer  (2)  Carbon  Layer  (3)  Heat 
Decomposition  or  Pyrolytic  Layer  (4)  Original  Materials  Layer 
(5)  Insulation  Layer 


Equations : 

3 

BY 


a 

37 


/,  ar\,  •  —  ar 

[k'wrm’'Cr'w =p' 


c,%-  X<J^X  +X 
ot 


ik'  BY  )  +  m*'  — AH,. 


(19) 


1>'W>  x  +  ^<7<JT  +X+X' 


In  this 


A'+X+A'/<7<A'c" 


(20) 


(21) 
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(22) 


P  ' c 


"  32  __  d  ( .  ,,,d T  \  X"<iY<.X>" 

’  a t  ~qy  \  a 7  )  0  0 


Original  conditions : 


T(Y,0)  =  g(Y)  =  Constant  Value 


Boundary  conditions : 

Carbon  layer  surface  Y  =  X 
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+  +  o  (23) 

,  07 

Common  boundary  surface  of  carbon  layer  and  thermal  decomposition 
or  pyrolytic  layer  Y  =  X  +  X 


•r=r'=r„ 


-k' 


bt'  _  _ , ar 

67  07 


(24) 


Common  boundary  surface  of  thermal  decomposition  or  pyrolytic 
layer  and  original  material  Y  =  X  +  X  +  X’ 


T  =  T"=Tn 
,dT' 


-k' 


07 


-k 


,,37" 

07 


(25) 


Common  boundary  surface  of  original  material  and  insulation 
material  Y  =  X' '0 


xf/ ~ x,n 

y^T"  _  0r"’ 

07  07 


(26) 


Insulation  material  back  surface  Y  =  X’ '  +  X'  1  ' 


or 


o 

67 

0T'1 ' 

h'  W~=E":'a  cr -r.‘D  (27) 
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3. 


EXPLANATION  OF  SYMBOLS 


An  explanation  of  the  different  meanings  of  symbols  is  seen 
in  Table  3 . 


TABLE  3  EXPLANATION  OF  SYMBOLS  WITH  DIFFERENT  MEANINGS 


CD  # 

X' 

X * 

*** 

X 

p" 

p" 

p’" 

\  i  ’Mm 

K' 

K" 

K 

<v 

V 

ct" 

Key:  Meaning  in  Decomposition  Zone  Algorithm  (4)  Original 
Material  Thickness  (5)  Insulation  Material  Thickness  (6)  Original 
Material  Density  (7)  Insulation  Material  Density  (8)  Original 
Material  Thermal  Conduction  Coefficient  (9)  Insulation  Material 
Thermal  Conduction  Coefficient  (10)  Original  Material  Specific 
Heat  (11)  Insulation  Material  Specific  Heat  (12)  Decomposition 
Zone  Material  Thickness  (13)  Original  Material  Thickness  (14) 
Insulation  Material  Thickness  (15)  Decomposition  Zone  Material 
Thickness  (16)  Original  Material  Density  (17)  Insulation  Material 
Density  (18)  Decomposition  Zone  Material  Thermal  Conduction 
Coefficient  (19)  Original  Material  Thermal  Conduction  Coefficient 
(20)  Insulation  Material  Thermal  Conduction  Coefficient  (21) 
Decomposition  Zone  Material  Specific  Heat  (22)  Original  Material 
Specific  Heat  (23)  Insulation  Material  Specific  Heat 


Explanations  for  symbols  with  the  same  meanings  are  seen  in 
Table  4. 
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TABLE  4  EXPLANATION  OF  SYMBOLS  WITH  THE  SAME  MEANINGS 


i)  ft  -§ 

0  ft  ft  £  X 

©  ft  .  ft 

Q  ft  ft  £  X 

X 

0  mmmm 

W 

mwAm&tfg)  _ 

X 

<d> 

e * 

mmmmm  & 

p 

s’" 

K 

e 

.  umatm 

cp 

<3)m  tm 

a 

Cp 

A  hc 

mc 

(2)  JfeJgfflSSlfc* 

Ahp 

(57)  "  " 

mt 

Tt 

j£fW.£S£t 

t 

qntt  fs]  . 

Tg 

#J8E3&K  ®> 

T 

T„ 

'  Qo  . 

<lR 

Key:  (1)  Symbol  (2)  Symbol  Meaning  (3)  Ablation  Layer  Thickness 
(4)  Carbon  Layer  Thickness  (5)  Carbon  Layer  Density  (6)  Carbon 
Layer  Thermal  Conduction  Coefficient  (7)  Carbon  Layer  Specific 
Heat  (8)  Thermal  Decomposition  or  Pyrolytic  Gas  Specific  Heat  (9) 
Carbon  Layer  Mass  Loss  Rate  (10)  Non-Carbon  Material  Loss  Rate 
(11)  Time  (12)  Temperature  (13)  Convective  Heating  Thermal  Flow 
Rate  (14)  Radiant  Heating  Thermal  Flow  Rate  (15)  Blocking  Effect 
Coefficient  (16)  Insulation  Material  Radiation  Coefficient  (17) 
Insulation  Material  Radiation  Coefficient  (18)  Material  Radiation 
Coefficient  (19)  Steven-Bortzmann  (phonetic)  Constant  (20) 
Combustion  Heat  Associated  with  Each  Unit  Mass  of  Carbon  (21) 

Heat  of  Thermal  Decomposition  or  Pyrolytic  Heat  (22)  Ambient 
Temperature  in  Cabin  (23)  Decomposition  Surface  Temperature  (24) 
Temperature  on  Common  Surface  of  Carbon  Layer  and  Decomposition 
Zone  (25)  Temperature  on  Common  Surface  of  Original  Material  and 
Decomposition  Zone 
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